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Abstract 

We derive closed expressions for the universal weak localization peak of the 

average conductance peak heights in Coulomb blockade quantum dots in the 

crossover from orthogonal to unitary symmetry. The scale for the crossover is 

independent of the number of channels in each lead, in contrast with the case 

of open dots. The functional form of the weak localization peak is independent 

of temperature. We also derive analytically the variance of the conductance 

peak heights as a function of the crossover parameter. 
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Quantum dots are formed in the interface between layers of semiconductors by the electro- 
static confinement of electrons in a two-dimensional electron gas. The transport properties 
(conductance) of such quantum dots can be measured by connecting them to external leads. 
In ballistic dots, the transport is dominated by scattering from the boundaries. Since dots 
with several hundred electrons are usually characterized by irregular shape, the scattering 
from the boundaries results in chaotic motion of the electrons. The conductance in such 
dots displays universal statistical fluctuations that are a signature of quantum chaos and 
that can be described by random matrix theory (RMT) [p]]. Earlier studies of quantum dots 
focused mostly on open dots, where the movement of the electron into the dot is allowed 
classically @. As the point contacts are pinched off, the dot becomes closed, i.e., effective 
tunnel barriers are formed in the interface between the dot and the leads. In such dots 
the conductance is dominated by resonant tunneling through electron resonances that are 
narrow compared with their average spacing. This leads to a series of narrow peaks in the 
conductance as a function of gate voltage, known as Coulomb blockade peaks. A statistical 
theory based on RMT for the fluctuations of the conductance peak heights was developed 
in Ref. || and was recently confirmed experimentally [[|,|5J . 

The statistics of the peak heights are universal but sensitive to the underlying symmetries. 
For conserved time-reversal symmetry (no magnetic field) the appropriate ensemble is the 
Gaussian orthogonal ensemble (GOE), while for broken time-reversal symmetry the statistics 
are described by the Gaussian unitary ensemble (GUE). The crossover from preserved to 
fully broken time-reversal symmetry is also universal. The distribution of the wavefunction 
intensity at a fixed spatial point was calculated by the supersymmetry method in Ref. 
||. However the fluctuations of the conductance in the crossover regime depend on the 
wavefunction fluctuations at several spatial points. This statistics and the corresponding 
distributions of the peak heights were derived as a function of a dimensionless crossover 
parameter in Ref. J7|. This crossover parameter is linear in the magnetic field that induces 
the symmetry breaking. 

One of the interesting phenomena associated with the onset of broken time-reversal sym- 
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metry in open dots is the suppression of weak localization. Weak localization in the absence 
of magnetic field originates from the constructive interference of pairs of time-reversed clas- 
sical trajectories, which enhances the return probability. This results in a smaller average 
conductance. When a magnetic field is applied, the constructive interference is destroyed 
and the weak localization is suppressed. In open dots, the semiclassical theory yields a 
Lorentzian shape for the weak localization peak ||. An exact expression for open dots was 
derived from RMT using the supersymmetry method |J, and approaches a Lorentzian for 
a large number of open channels. 

Weak localization also occurs in almost closed dots, where the average conductance peak 
height was found to be smaller in the GOE case than in the GUE case H. On the other 



hand, no weak localization effect is found in the average of the conductance minima |[T0| . pT 
In this paper we derive closed expressions for the weak localization peak of the conductance 
maxima in closed dots as a function of the crossover parameter and for any number of 
equivalent channels in the leads. We also study the temperature dependence of the weak 
localization peak and find that, when measured in units of the average GUE conductance, 
the peak is temperature-independent (assuming no phase breaking). We find that the width 
of the peak is only weakly dependent on the number of channels in the leads, in contrast 
to the case of open dots, where the width of the peak is proportional to the square root 
of the number of open channels. We also derive closed expressions for the variance of the 
conductance peak as a function of the crossover parameter. 

For almost closed dots, the resonances in the dots are isolated and their average width 
is small compared with the mean level spacing A. At temperatures that are low compared 
with the mean level spacing (kT <C A), only one resonance contributes to a given Coulomb 
blockade peak, and the conductance peak amplitude G is G = where ]T2] 

i r*r r 

Here T 1 ^ is the width of a resonance to decay into the left (right) lead and F is the average 
width per channel. Note that the quantity g is dimensionless and temperature- independent. 



In general we assume that the left (right) lead has A l( - r ^ open channels such that T 1 ^ = 
J2 C ITc^I 2 ) where 7^ is the partial amplitude to decay into channel c on the left (right). In 
i?-matrix theory, the partial width amplitude 7 C can be expressed as a scalar product (defined 
over the dot-lead interface) of the resonance wavefunction ip and the channel wavefunction 
4> c through 7 C = (4> c \ip) oc fdS</>*(r)if)(r). We assume for simplicity that the A channels are 



uncorrelated and equivalent, i.e. 7*7 C ' = T6 cc i. All average partial widths are assumed to 
be independent of the magnetic field. We also assume that the dot has symmetric leads, so 
that the average partial widths are the same in the left and right leads. 

For a ballistic dot with chaotic dynamics of the electrons, we model the Hamiltonian H 
of the dot by an ensemble of random matrices. When we apply a magnetic field to the dot, 
the ensemble describing the crossover from GOE to GUE is |13|JT^ 



H = S + iaA, (2) 

where S and A are uncorrelated symmetric and antisymmetric real matrices, respectively, 
chosen from Gaussian ensembles of the same variance, and a is a real parameter. The 
transition parameter describing the crossover is given by the ratio of the root-mean-square 



(rms) of a typical symmetry-breaking matrix element to the mean-level spacing \T5\ A = 
(H 2 break) 1 ^ 2 / A = ocyN / V . Alternatively 2ir\ = y th j T m { x , where th = h/A is the Heisen- 



berg time and r mix is the mixing time defined by the spreading width h/r mix = 27iH 2 break / A 
due to the interaction that breaks time-reversal symmetry . The spectral properties of the 
ensemble (|2|) make the complete crossover for A ~ 1. The average dimensionless conductance 
peak height g(\) can be calculated by averaging (HD over the ensemble (H). 

The statistics of the partial width amplitudes 7 were derived in Ref. 0. The components 
ij^n of the resonance wavefunction (in a fixed basis) are decomposed into their real and 
imaginary parts ip^ = ip^R + iip^j in a principal frame in the complex plane (determined 
by J2^=i'4 ) ^R' t P^i = 0), and a parameter t is defined by t 2 = £^ tp 2 ^ ^ i^R- At a fixed 
value of t, the real and imaginary parts of a finite number A of components (A N) are 
independent Gaussian variables. Except in the GOE and GUE limits (where t — and 



t = 1, respectively), t is not sharp but fluctuates according to a known distribution []6|,|7|,[16| 



Px(t) = ^l^AV^-V*) 2 L(\) + 



1 / 1\ 2 1 
t + 



4 V ' tj 2vr 2 A 2 



[1 - 0i(A)] , (3) 



where X (A) = / e^^^dy. These fluctuat ions of t are responsible for the long-range 
o 



correlations of wavefunction intensities in the crossover from GOE to GUE [17 . 

Applying a similar decomposition for the partial width amplitudes 7 C = 7 c r + i^ c i = 

(4> c \ipn) + i(<fi c \il>i) in the principal frame of ip, it is found that at fixed t, 7# and 7/ 

are Gaussian statistically independent variables with 7^ = t 2r y 2 R . The joint distribution of 

the partial widths is obtained by averaging the conditional distribution P("y\t) over the t 
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distribution, ^(7) = / P x (t)P(j\t)dt. 



Using Eq. (|l]) and the statistical independence of the partial widths in the left and right 
leads at fixed t (i.e., P(T l ,T r \t) = P(T l \t)P(T r \t) where P(. . . \t) are conditional probabili- 
ties), the average conductance peak height can be written as 

^(A) = ^°°rfrQ( 1 )(r,t)/r^ (4a) 

Q (1) (I\ t) = f r dT l T l P(T l \t)(T - T l )P(T - T l \t) , (4b) 
J o 

where (...) denotes here and in the following an average over the t distribution P\(t) given 
by Eq. (^|). The function Q^(T,t) is the convolution of T l P(T l ) with itself, and its Laplace 
transform is Q^\s, t) = (dP/ds) 2 where P(s, t) = dTe~ sT P(T\t) is the Laplace transform 
of the width distribution. P(s, t) is easily calculated using the Gaussian nature of the partial 
width amplitudes at fixed t and T = J2c (|7cr| 2 + |7e/| 2 )- Assuming one-channel symmetric 
leads (i.e. T l = T l = T) and measuring T in units of T, we find Pi(s, t) — (1 + s 2 /x 2 + 2s) 
where x = (t _1 +t)/2, and thus Qi\s,x) = x 2 (s + x 2 ) 2 (s 2 + 2sx + x 2 )~ 3 . The inverse Laplace 
transform Qi(T, t) is calculated by the residue theorem using the poles (of degree three) of 
Qi(s, x) at s — t and s = 1/t. After performing the integral in Eq. (^a[), we find the average 
conductance as a function of A to be 

^ + ((r^)m- + i)V 



The GOE and GUE limits are obtained for t — > and t — > 1, respectively, where <?goe = 1/4 
and (fcUE =1/3 0. 

The above method can be generalized for any number A of equivalent and uncorre- 
cted channels (in each of the symmetric leads). Measuring T in units of the average 
resonance width per channel, we find P^(s,t) = (1 + s 2 /x 2 + 2s) A ^ 2 and Q^\s,t) = 
A 2 x 2A (s + x 2 ) 2 (s 2 + 2x 2 s + x 2 ) _A_2 . Using (§) we find that #a changes from A 2 /2(A + 1) in 
the GOE limit (x — > oo) to A 2 / (2A + 1) in the GUE limit (x — > 1). In the crossover regime 
we find that <?a(A) has the form 

( 2A N 



^-2(aTi) + H a+1 a 



V A-I y 



2A 



2t 2 4=J A-t 2 \ 2n 



In t - a n 



1 1 n=0 V 1 



• (6) 



The coefficients a n in Eq. (^j) can be determined from the known GUE limit (t — > 1) of g~\. 
For the limit on the r.h.s. of (0) to exist when t — > 1, the coefficients a n for n = 0, 1, . . . , A — 1 
must be equal to the corresponding expansion coefficients of the function 2t 2 / (1 — t A ) lnt in 
powers of [(1 — t 2 )/t] 2 , and are therefore independent of A. The value of g\ in the GUE limit 
also determines the coefficient a\ to be a\ = (— ) A+1 (A!) 2 /2((2A + 1)!. Since the coefficients 
independent of A, we conclude 

fn!) 2 



\n+l 



(7) 



n v ' 2(2n + l)! ' 

The function 5^a(A) = g^ VE — <7a(A), which is peaked at A = (GOE limit) and ap- 
proaches zero for A — > ±oo, is known as the weak localization peak. Its analytic form is 
plotted in Fig. [I] (in units of the average GUE conductance peak) for A = 1, 2, 3 and 5. 
The size of the peak is 5g A (0) = ^ UE - ^ OE = A 2 /[2(A + 1)(2A + 1)]. For open dots 
with A equivalent channels (with transmission coefficient 1) in each lead the size of the 
weak localization peak is A/[2(2A + 1)] ||, so that the weak localization peak in closed 
dots (measured in units of ^ JrL) is smaller by a factor of A/(A + 1) than the one in open 
dots (measured in units of e 2 /h). In units of the average GUE conductance height, we find 
^a(0)/^a UE = 1/[2(A + 1)] for closed dots, compared with 1/(2A + 1) for open dots. 
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It is interesting to compare the dependence of the width of the weak localization peak 
on the number of channels in open versus closed dots. The right inset of Fig. [I] shows 
the full width at half maximum (FWHM) of the weak localization peak as a function of 
the number of channels for both closed (solid circles) and open (open circles) dots. The 
FWHM in closed dots is approximately independent of the number of channels, unlike the 
case of open dots where the FWHM increases as y/X (for large A). In open dots, the 
crossover scale in the conductance is determined by the competition between the mixing 
time Ti/r m i X = 27rAA 2 (associated with the breaking of time-reversal symmetry) and the 
decay time of resonances in the dot h/Td ec = (A/27r)A. The crossover in the conductance 
occurs when T dec /T mix = An 2 X 2 /A ~ 1 i.e., for X^ en ~ \ZA/2ty. The decay time in an open 
dot is shorter than the Heisenberg time by a factor of A, and for one open channel these two 
times are equal. For closed dots on the other hand, the decay time is much longer than the 
Heisenberg time (since the average level width is much smaller than A), irrespective of the 
number of channels. In such a case the relevant time scale that competes with the mixing 
time is the Heisenberg time, and the crossover occurs on a scale of \f° sed ~ 1, independent 
of the number of channels and comparable to the value of A cr for an open dot with A = 1. 
This is confirmed by the exact RMT results of Fig. where the full width at half maximum 
(FWHM) of the weak localization peak is AA ~ 0.4, irrespective of the number of channels 
in each lead. 

If time-reversal symmetry is broken by a magnetic field, then A = $/$ cr , where $ cr 
is the crossover flux (defined as the flux where A = 1). It is possible to estimate $ cr 
semiclassically assuming single-particle chaotic motion PJl8fl. The rms of the electron action 
difference between a pair of time-reverse orbits is given by [(0 2 ) 1 / 2 /^4]($/$o) (in units of 
h), where (O 2 ) 1 / 2 is the rms area accumulated by the electron and A is the area of the 
dot. Since the area accumulation in a chaotic dot is diffusive, we have for an open dot 
(O 2 ) 1 / 2 1 A = ceg 1 / 2 ' \J T cr / Tdec, where r cr is the time scale for the electron to cross the dot, 
and a g is a geometrical factor that is dot specific (Tj|. For the breaking of time reversal 
symmetry, we require the rms action difference to be of order 1, leading to $°^ en /$o = 
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na g 1 ^ 2 y/h.\Jr cr /rH- For the closed dot, on the other hand, th replaces Td ec to give a A- 
independent crossover flux of §f° sed /&o = ^Qg^ 1 ^ 2 \/ T cr/ T H- A weak localization of the 



average conductance peak height was observed experimentally in semi-open dots |TT| with 
an FWHM of ~ 6.2 mT, corresponding to B cr « 14.9 mT or $ cr « 2<3>o. The desymmetrized 



stadium billiard in a uniform magnetic field gives $ cr ~ |$o HU with about the same number 
of electrons (~ 1000). Thus the experimental result is about three times larger than the 
billiard model estimate, similar to the discrepancy estimated for the correlation field J?|. 
The exact discrepancy is uncertain because of the unknown geometrical factor of the dot 
used in the experiment. The discrepancy does not indicate a breaking of RMT (which only 



predicts the universal form of the correlator [19| and not its magnetic flux scale), but rather 



that effects beyond single-particle dynamics are important. Indeed, a model that takes into 
account electron-electron interactions in the dot gives a correlator whose form is consistent 
with the RMT correlator, but with a correlation field that is larger by about a factor of 
three as compared with a non-interacting model ^0] . 

It is also possible to calculate in closed form the variance of the conductance as a function 
of A. The second moment of the conductance is given by g 2 (X) = ^/ °° dTQ^(T, t)/T 2 ^, 
where is the convolution of (T l ) 2 P(T l \t) with itself. The Laplace transform is then 
calculated from Q {2 \s,t) = [d 2 P{s\t)/ds 2 ] . For one-channel leads we find 



lnt+1 ' 1 ' 



1 - 1 2 12 V t J 27 V t 

g 2 increases from 3/16 (GOE) to 1/5 (GUE), while the variance <J 2 (g) = g 2 — g 2 decreases 
from 1/8 to 4/45. Closed expressions for g\(X) can be similarly derived for any number of 
equivalent channels A. The left inset of Fig. [I] shows the analytically calculated a(g)/g as a 
function of A for A = 1,2,3 and 5 channels. We see that the distribution becomes sharper 
in the crossover from conserved to broken time-reversal symmetry. In the GOE and GUE 
limits we have 

gte)/ « = < (CUE) " ( ) 

2A 2 (2A+3) J l^urvj 



(8) 



For large A, this quantity decreases in the crossover from 1/vA (GOE) to l/v2A (GUE), 
i.e., effectively doubling the number of channels. 

At finite temperature, several levels i contribute to the conductance peak with known 



thermal weights Wi(T) [T2|]2]J which are determined by the single-particle spectrum and the 
number of electrons in the dot. The dimensionless peak height is given by g = J2iWi(T)gi, 
where Qi are the level conductances as in Eq. ([[]), but with r' ,r depending now on the level 
i. In calculating the statistics of the peak heights at finite T, a good approximation is to 
neglect the fluctuations in the dot's spectrum, and Wi(T) are then fixed. Since g«(A) are 
independent of i, the weak localization peak factorizes 5g(X,T) « Wi(T)} Sg(\), where 
Sg(X) is the weak localization peak at temperatures much smaller than A. In particular, 
since g E (T) « [Y^,i w i{T)} g GlJE , the scaled weak localization peak 5g/g GlJE is temperature- 
independent (assuming no phase breaking). 

Another quantity of interest is the size of the weak localization peak relative to the 
rms of the (GUE) peak-heights distribution. For dots with one-channel leads, we find 
5g(0)/a GlJE (g) « (y/E/8) [(J2 Wi) 2 / W V^^ 2 which increases with temperature. At high 
temperatures, the number of levels that contribute to a given peak height is of the order 
of T/A, and 5g(0)/aG\jE(g) oc (T/A) 1 / 2 . At low temperatures it is difficult to measure the 
weak localization peak because its size is comparable to the conductance fluctuations. In 
open dots, the large number of open channels enhances the peak relative to the fluctuations 
by y/X. In closed dots (with A = 1), the temperature can be used instead to enhance the 
peak relative to the background fluctuations. 

In conclusion, we have derived in closed form the weak localization peak of Coulomb 
blockade peak heights in quantum dots. In the absence of phase breaking, this peak is 
temperature-independent up to an overall scale. 

We acknowledge CM. Marcus for helpful discussions. This research was supported in 
part by DOE Grant DE-FG02-91ER40608. 
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FIGURES 
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FIG. 1. The weak localization peak in closed dots (measured in units of the average GUE 
conductance) Sg(X)/g GVE versus the crossover parameter A. Shown are the cases with A = 1,2,3 
and 5 channels (peak is decreasing as A increases). Right inset: The FWHM of the weak localization 
peak AA as a function of the number of channels A in each lead for closed and open dots (solid 
and open circles, respectively). Left inset: The ratio a(g)/g between the rms and average values of 
the conductance peak as a function of the crossover parameter A for Coulomb blockade dots with 
A = 1, 2, 3 and 5 channels (a(g)/g is smaller for larger A). 
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